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Abstract 



The four LEP collaborations, ALEPH, DELPHI, L3 and OPAL, have searched for pair-produced 
charged Higgs bosons in the framework of Two Higgs Doublet Models (2HDMs). The data of the 
four experiments are statistically combined. The results are interpreted within the 2HDM for 
Type I and Type II benchmark scenarios. No statistically significant excess has been observed 
when compared to the Standard Model background prediction, and the combined LEP data 
exclude large regions of the model parameter space. Charged Higgs bosons with mass below 
80 GeV/c 2 (Type II scenario) or 72.5 GeV/c 2 (Type I scenario, for pseudo-scalar masses above 
12 GeV/c 2 ) are excluded at the 95% confidence level. 
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1 Introduction 



More than twelve years after the end of data-taking at LEP, it is still useful to add to the LEP 
legacy the outcome of the searches for a charged Higgs boson. In fact, the charged Higgs boson 
searches at hadron colliders are restricted to specific decay channels, thus having (yet) some 
difficulties to provide limits in the charged Higgs boson mass independent of tan/3 (the ratio 
of the two Higgs doublet vacuum expectation values). 

Since the previous communication by the LEP working group for Higgs boson searches 
(LEPHWG) on charged Higgs boson searches, in 2001 PQ, the LEP experiments have published 
their final results on these searches and have, in some cases, also added searches for new final 
states not previously considered. The four LEP collaborations have searched for charged Higgs 
bosons in the framework of Two Higgs Doublet Models (2HDMs). Based on the final results 
obtained by ALEPH 0, DELPHI [3,4], L3 and OPAL [6,7], the LEPHWG has performed 
a statistical combination of the data taken at centre-of-mass energies, a/s, from 183 GeV to 
209 GeV. The total luminosity used in this combination is 2.6 fb~ . 

The existence of a pair of charged Higgs bosons is predicted by several extensions of the 
Standard Model (SM). In 2HDMs, the charged Higgs couplings to the photon and the Z boson 
are completely specified in terms of the electric charge and the weak mixing angle, 8\y, and 
therefore, at tree level, the production cross-section depends only on the charged Higgs boson 
mass. Higgs bosons couple proportionally to the particle mass and therefore decay preferentially 
to heavy particles, but the precise branching ratios may vary significantly depending on the 
model. Two scenarios are considered in this paper. The first one effectively allows the charged 
Higgs boson to decay to fermions only, which is the case in type II 2HDM [9] for not too 
small values of (the neutral CP-odd A boson mass) or tan/3. In this model the isospin +| 
fermion-couplings to the charged Higgs boson are proportional to 1/ tan /3, while the isospin — | 
fermion-couplings are proportional to tan (3. This scenario is treated in Section 2. In the second 
scenario, type I 2HDM [10], all fermions couple proportionately to 1/ tan/3. Consequently, the 
second scenario effectively allows the charged Higgs boson to also decay into gauge (possibly 
off-shell) bosons and Higgs bosons (see Section 3). 

At LEP, charged Higgs bosons are pair-produced, via s-channel exchange of a Z° boson (real 
or virtual). Electroweak precision measurements set indirect bounds on the mass of the charged 
Higgs boson regardless of its decay branching ratios. The difference between the measured decay 
width of the Z° (rz) and the prediction from the SM sets a limit on any non-standard (non 
SM) contribution to Z° decay. The Z° decay width has been measured precisely during the 
first phase of LEP (LEP-1). The final LEP result [UJ set the limit r nonSM < 2.9 MeV/c 2 (95% 
C.L.), which translates to mn± > 39.6 GeV/c 2 (95% C.L.). Direct searches during the LEP-1 
period set a lower bound for the charged Higgs boson mass at 44.1 GeV/c 2 at 95% C.L. for 
type II 2HDM [12-15]. The combination in this paper is performed for charged Higgs boson 
masses of 43 GeV/c 2 or larger, since the region below 43 GeV/c 2 has been covered by individual 
experiments. 

For this combination of data, the cross-sections (and branching ratios for type II 2HDM) 
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are calculated within the HZHA program package [TB], and the branching ratios of the charged 
Higgs boson in type I 2HDM are taken from Ref. [T7] . 

The input from the four experiments [2-7] which is used in the combination procedure is 
provided on a channel-by-channel basis. The word "channel" designates any subset of the data 
where a Higgs boson search has been carried out. These subsets may correspond to specific 
final-state topologies, to data sets collected at different centre-of-mass energies or to the subsets 
of data collected by different experiments. Table [T] shows a summary of all channels available 
for this combination. It amounts to 22 channels from ALEPH, 43 from DELPHI, 12 from L3 
and 45 from OPAL. 

Each experiment generated and simulated the detailed detector response in Monte Carlo 
event samples for the Higgs signal and the various background processes, at centre-of-mass 
energies of 183, 189, 192, 196, 200, 202, 204, 206, 208 and 209 GeV to estimate background 
and signal contributions in the data collected between 1997 and 2000. Particular care has been 
taken when simulating the four-fermion background, especially from W-pair background, using 
the most advanced codes available at that time. ALEPH used KORALW [18] as the generator 
and RACOONWW 02] and YFSWW [20] for the cross-section calculation, while DELPHI used 
WPHACT [21 j, L3 YFSWW and OPAL GRC4F [22] and KORALW. Other generators were 
used for systematic studies. Furthermore, each of the four experiments used different values 
for the W mass in these background simulations (respectively 80.45, 80.40, 80.356 and 80.33 
GeV/c 2 for ALEPH, DELPHI, L3 and OPAL), while the LEP combined measured value is 
80.376 GeV/c 2 [25] . thus introducing an additionnal source of systematic uncertainty. 
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Table 1: Overview of the searches for charged Higgs bosons performed by the four LEP experiments, 
whose results are used in this combination. Where relevant, tha varies from 2mf, to "ijj± • 
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The statistical procedure adopted for the combination of the data and the precise definitions 
of the confidence levels CL S , CL s+ t, by which the search results are expressed, have been 
described previously [21]. The main sources of systematic uncertainty affecting the signal 
and background predictions are included, using an extension of the method of Cousins and 
Highland [25] where the p- values are averaged over a large ensemble of Monte Carlo experiments. 
The correlations between search channels, LEP collision energies and individual experiments 
have not been taken into account, but these correlations are estimated to have only small effects, 
about 500 MeV/c 2 , to the final results. 



2 Combined searches in the framework of type II 2HDM 

In 2HDMs there are five physical Higgs bosons: the CP-even h and H, the CP-odd A and the 
charged Higgs bosons, H 1 * 1 . The various types of the 2HDMs differ by their couplings to the SM 
fermions. In type II 2HDM [9], one Higgs doublet couples to up- type fermions and the other 
to down-type fermions. The Higgs sector of the Minimal Supersymmetric Standard Model 
(MSSM) is a particular case of such models. In the MSSM, at tree-level, the H ± is constrained 
to be heavier than the W boson and the radiative corrections to the charged Higgs mass are 
positive, except for very specific parameter choices. Thus, experimentally finding evidence of 
a charged Higgs boson with mass below the W boson mass would set very strong constraints 
on the MSSM parameters. However, in the following we will concentrate on the general type II 
2HDM without any supersymmetric assumptions. Results on the search for neutral MSSM 
Higgs bosons can be found in [26J. 

For the charged Higgs masses accessible at LEP energies, the decays into r + z/ r and cs (and 
their charge conjugates) are expected to dominate. The searches are carried out under the 
assumption that the two decays H + — >cs and H + — >t + u exhaust the H + decay width, but the 
relative branching ratio is free. This assumption is valid as long as wa is larger than 60 GeV/c 2 
(MSSM case) or tan/3 is large. Thus, the searches encompass the following H + H~ final states: 
(cs)(cs), (r + 'v)(t~ v) and the mixed mode (cs)(t~u) or (cs) (r + z/). The combined search results 
are presented as a function of the branching ratio Br(H + — >t + u). 

Details of the searches done by the individual experiments can be found in Refs. [2-7]. Two 
features in these analyses are worth noting: the main background is W pair production, which 
is partly irreducible, and the reconstructed mass is one of the discriminant variables used in 
the final hypothesis testing in the two channels where this is relevant (mixed and hadronic 
channels). The results from the four LEP experiments are summarised in Table [21 together 
with the 95% C.L. observed and median expected lower limits on the charged Higgs boson 
mass. The mass limits are quoted separately for Br(H + — >r + v) = 0, 1, and independently of 
the charged Higgs decay. 

In a first step, the statistical combination software was run separately on the data provided 
by the four collaborations and the results (shown in Figured]) compared to the published results 
of each. The differences between this check and the published results reflect the differences 
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Table 2: Individual search results for the e + e — t-H + H fermionic final states. All limits are given 
at the 95% C.L.. (*) The OPAL selection is mass- dependent; the numbers given here are for mjj± = 
80 GeV/c 2 . 

between the statistical methods used by the four collaborations. The biggest difference has 
been found for the expected limit from OPAL at Br(H + — >r + v)= 1 and amounts to 600 MeV/c 2 . 
This difference is compatible in size with the estimated effect of about 500 MeV/c 2 of not taking 
into account the correlations between systematic uncertainties. All mass limits have thus been 
rounded down to the nearest half a GeV/c 2 . 

Combining the results from the four experiments, a scan in the branching ratio Br(H + — >t + u) 
versus charged Higgs boson mass plane has been performed, and the limit-setting procedure 
was repeated for each scan point. This two-dimensional scan was performed with the following 
ranges and steps: wh± from 43 to 95 GeV/c 2 with 1 GeV/c 2 steps, and Br(H + — >t + u) from 
to 1 with 0.05 steps. 

Figure [2] shows the observed background confidence level CLb as a function of mn± and 
Br(H + — >r + u). The observed confidence level is everywhere within ±2cr of the background 
prediction, except for three small regions, as shown in Figure 2. Such regions result from the 
combination of small excesses, as compared to the background expectation, observed by two or 
three experiments. The first two mass regions, around 43 and 55 GeV/c 2 , are due to a slight 
excess of data in the hadronic channel, while the third one, around 90 GeV/c 2 , arises from an 
excess of events in the mixed channel. Table [3] gives the combined CLb together with the values 
from each experiment for these three domains (the values chosen for Br(H + — »r + z/) are given in 
the second column). 

The combined results for the Type II 2HDM are summarised in Figure [3l which shows the 
expected median and observed mass limits, while the contribution of each of the three decay 
channels to the overall limit is presented in Figure HI It is worth noting that: 

• the purely leptonic channel alone excludes charged Higgs masses above the W mass, down 
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Table 3: Combined and individual CL b values for the three mass points with a deviation from expec- 
tation larger than 2 a. All values, obtained with the statistical procedure of the overall combination, 
compare well to those published by the experiments. (*) ALEPH and L3 did not provide inputs for 
this mass. 

to Br(H + — >t + u) around 0.45. In this channel, the mass of the Higgs boson cannot be 
reconstructed, due to the presence of two neutrinos in the final state. As a consequence, 
the W boson pair background is diluted and the analysis is sensitive up to y/s/2. The 
limit drops rapidly for Br(H + — >t + u) below 0.45, due to a rapid decrease of the signal rate 
in this final state; 

• the mixed channel alone cannot exclude charged Higgs mass values up to the W mass, 
even when it contributes maximally, for Br(H + — >t + v) = 0.5. For this value, the observed 
limit is only slightly above 79 GeV/c 2 , due to the large e + e~— )-W + W~ background. This 
channel has the best coverage in terms of Br(H + — >t + u), as shown in Figure 4; 

• the hadronic channel is the most difficult one; for masses close to the W mass, the 
sensitivity is reduced due to the large e + e~— >-W + W~ background. The sensitivity at 
higher masses is improved (a gain of 10 GeV/c 2 on the expected limit), and the observed 
limit as well (note the excluded "island" at Br(H + — >t + u) close to zero) with respect to 
the results of individual experiments; 

• the difference between the expected and observed limit seen in Figure [3]for Br(H + — >t + v) 
from 0.35 to 0.85 results from the excess of observed events already mentioned (see Figure 
2) in the mixed channel above m H ± = 84 GeV/c 2 . 

The combined 95% C.L. mjj± lower limits are listed in Table H] for Br(H + — >t + u)=0, 1, 
together with the limit that is independent of the fermionic decay mode. Taking the lowest of 
the observed limits from Table HI we quote a 95% C.L. lower bound of 80 GeV/c 2 for the mass 
of the charged Higgs boson in type II 2HDM under the assumption of pure fermionic decays 
of the charged Higgs boson. Thus the hypothesis of a charged Higgs boson degenerate in mass 
with the W boson is not excluded at the 95% confidence level with LEP data. The limits 
around the W mass are very sensitive to the modelling of the W pairs background. Taking the 
uncertainties in the background modelling into account results in a downward shift of these 
limits by 600 and 500 MeV/c 2 for Br(H + — >t + v)=0 and 0.5, respectively. 

Figure shows the 95% C.L. upper bound on the e + e~— >-H + H~ cross-section (with ± la 
and ± 2a bands) for four values of Br(H + — >r + u), namely 1, 0.5, 0.2 (which corresponds to 
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H ± mass limit (95% C.L.) 
in GeV/c 2 


Br(H+->r+z/)=0 
Expected limit (median) 
Observed limit 


88 
80.5 (*) 


Br(H+-^r+z/) = l 
Expected limit (median) 
Observed limit 


93.5 
94 


Any Br(H+^r+z/) 
Expected limit (median) 
Observed limit 


79.5 
80 



Table 4: The combined 95% C.L. lower bounds on the mass of the charged Higgs boson, expected 
and observed, for fixed values of the branching ratio Br(R + — >t + v) and for any Br(H + ^-T + u). All 
mass limits have been rounded down to the nearest half a GeV/c 2 to take into account the effect of 
neglecting the correlations between systematic uncertainties. (*) The interval from 83 to 88 GeV/c 2 
is also excluded at the 95% C.L.. 

the weakest limit) and 0. The thick black curve is the 2HDM tree-level prediction for that 
cross-section. 
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Figure 1: Type II 2HDM: the 95% C.L. bounds on mjj± as a function of the branching ratio 
Br(Yi + ^T + u), for each of the four LEP experiments separately. The expected exclusion limits are 
indicated by the thin solid line and the observed limits by the thick solid line. The shaded area is 
excluded at the 95% C.L. or higher. 



7 




50 60 70 80 90 



Charged Higgs mass (GeV/c) 

Figure 2: Type II 2HDM: contours based on the observed p-values CLf, as a function of mjj± and 
the branching ratio Br(H + — >t + u), indicating the statistical significance, N a , of local departures from 
the background expectation. The black solid line indicates the change of sign of this significance, i.e. 
where there is a transition from excess to deficit. 
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Figure 3: Type II 2HDM: excluded regions in the Br(H + ^-r + u) vs m H ± plane, based on the combined 
data collected by the four LEP experiments at centre- of-mass energies from 183 to 209 GeV. The 
shaded area is excluded at the 95% or higher C.L.. The expected exclusion limit (at the 95% C.L.) is 
indicated by the thin solid line and the thick dotted line inside the shaded area is the observed limit at 
the 99.7% C.L.. 
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Figure 4: Type II 2HDM: regions in the Br(H + ^t + u) vs mjj± plane excluded at the 95% or higher 
C.L., based on the combined data collected by the four LEP experiments at centre- of-mass energies 
from 183 to 209 GeV, for each of the three decay channels separately. The solid (dashed) lines are the 
observed (expected) limits. 
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Figure 5: Type II 2HDM: the 95% C.L. upper limits on the production cross-section as a function 
°/ m H± for four different values of the branching ratio Br(a + ^T + u), combining the data collected by 
the four LEP experiments at centre- of-mass energies from 183 to 209 GeV. The solid lines represent 
the oserved exclusion limits, while the expected exclusion limits are indicated by the dashed lines. The 
shaded bands represent the =L1<7 and ±2<r excursions around the expected limits. The intersections 
of the curves (solid or dashed) with the thick line showing the theoretical (tree-level) charged Higgs 
cross-section represent the (observed or expected) 95% C.L. lower limits on the charged Higgs boson 
mass. 



10 



3 Combined searches in 2HDM of type I 



An alternative set of models, type I 2HDMs [ID] , assume that all fermions couple to the same 
Higgs doublet. In this case all fermions couple proportionately to 1/tan (3 to the charged Higgs 
boson and fermionic decays are suppressed for medium to large tan (3 values. Consequently, if 
a neutral Higgs boson $ (representing either A or the lightest CP-even scalar h) is sufficiently 
light, the decay to W*$ can be dominant even in the range of charged Higgs masses of interest 
at LEP (where W* indicates an off-shell W boson). While searches for a CP-even neutral Higgs 
boson exclude such a particle for masses below 82 GeV/c 2 independently of its decay [27], the 
existence of a light CP-odd neutral Higgs boson, A, is not excluded by experiment [28]. Hence, 
the search for the process H 1 * 1 — >W*A is fully justified. Figure shows the predicted branching 
ratios of the charged Higgs bosons for various choices of parameters of type I models. For all 
kinematically allowed values of the A mass, tua, the possible charged Higgs boson decays are 
predominantly fermionic for low tan (3 and predominantly bosonic for high tan /3. Between these 
two extreme cases, the branching ratios change rapidly as a function of tan (3 (between typically 
0.1 and 10) and slower as a function of m^, appearing earlier in tan/3 for lower m^- The ratio 
between the two competing fermionic decays irv over cs) is almost independent of the charged 
Higgs boson mass (see lower part of the figure), as expected from the Yukawa coupling which 
only depends upon the masses involved. 

To cover the possibility of a light A boson the final states W*AW*A and W*Ar~z7 r were also 
searched for by DELPHI [3 J and OPAL [6]. The channel W*Acs was not considered because 
its contribution is expected to be small for all tan (3. The A boson was searched for through its 
decay into two b-jets, restricting the A mass to be above 12 GeV/ c 2 . Type I models are explored 
through the combination of all five decay channels, namely the final states cscs, csrz/, t + vt~v, 
W*AW*A and W*Ar - z7 T (and their charge conjugates). The combination of the experimental 
search results is performed for branching ratio values predicted by the model as a function 
of tan/3 and m^. Where there was a possible overlap between two search channels, the one 
providing less expected sensitivity was ignored to avoid double counting. This is the case in the 
intermediate region in tan f3 for purely hadronic channels (W*AW*A and cscs) on the one hand 
and the semi-leptonic channels (W*Ar~z7 r and cst~P) on the other. A three-dimensional scan 
was performed with the following ranges and steps: mn± from 43 to 95 GeV/c 2 in 1 GeV/c 2 
steps, m\ covering 12 GeV/c 2 , then 15 to 75 GeV/c 2 in 5 GeV/c 2 steps, and tan/3 from 0.1 to 
100 in steps of 0.2 in log(tan/3). 

Figure [7J shows the observed CLj,, for four values of and two values of tan/3. A slight 
excess for low and intermediate A masses in the high tan/3 region where the bosonic decays 
dominate is observed, resulting in observed limits generally weaker than expected (see Figure [8]). 
Three main features are visible in Figure [HJ two plateaux and a valley between them: 

• the first plateau, at low tan/3, corresponds to the case when the fermionic channels 
dominate. Both expected and observed limits are above 86 GeV/c 2 ; 

• the valley is somewhat of an artefact. It is due to the conservative approach of considering 
only the most sensitive channel when two overlapping channels contribute. The difference 
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between expected and observed mass limits reaches 4.5 GeV/c 2 in the extreme case (when 
tan/3 = 1.6 and m A = 12 GeV/c 2 ); 

• the second plateau, at high tan /3, corresponds to the case when the bosonic channels 
dominate. The small excess seen in Figure [7] corresponds to a small difference between 
expected and observed charged Higgs mass limits, which is always less than 2.2 GeV/c 2 . 

Table [5] summarizes these results. For low tan/3 (tan/3 below 0.5) where the bosonic contri- 
bution is vanishingly small, the ihh± lower limits (above 86 GeV/c 2 ) are almost independent of 
mA- On the other hand, for high tan/3 (equal to or greater than 10) where the bosonic chan- 
nels dominate, the sensitivity is maximal for intermediate A masses (mA around 50 GeV/c 2 ). 
Outside the valley, the limit is always above 84 GeV/c 2 . Finally, the lowest limits always corre- 
spond to the cases in the valley, thus depending both on tan/3 and m^. The lowest (observed) 
limit is 72.5 GeV/c 2 , for tan/3 = 1.6 and m A = 12 GeV/c 2 . This limit rises to 76.5 GeV/c 2 
for rriA = 20 GeV/c 2 and the difference between expectation and observation is reduced to 1 
GeV/c 2 . 



m A 


tan/3 = 0.1 


tan (3= 1 


tan/3= 10 


tan (3= 100 


minimum 


12 


86.0 (86.0) 


73.5 (77.0) 


83.5 (86.0) 


84.0 (86.0) 


72.5 (77.0) 


20 


86.5 (86.0) 


76.5 (77.5) 


85.5 (87.0) 


85.5 (87.0) 


76.5 (77.5) 


30 


86.5 (86.5) 


80.0 (79.5) 


87.5 (89.0) 


87.5 (89.0) 


78.0 (79.5) 


50 


86.5 (86.5) 


84.0 (84.0) 


89.0 (90.0) 


89.5 (91.0) 


81.0 (80.5) 


70 


86.5 (86.5) 


86.5 (86.5) 


83.5 (83.5) 


89.0 (90.5) 


81.0 (81.0) 


minimum 


86.0 (86.0) 


73.5 (77.0) 


81.5 (81.0) 


81.0 (81.0) 


72.5 (77.0) 



Table 5: Observed lower limits on the charged Higgs mass in GeV/c 2 at 95% C.L. for different values 
°f m A (in GeV/c 2 ) and tan/3. The expected median limits are shown in parentheses. The last column 
(last row) show the weakest limit for a fixed A mass and any tan (3 (for a fixed tan (3 and any A mass). 
The mass limits have been rounded down to the nearest half a GeV /c 2 to take into account the effect 
of neglecting correlations between systematic uncertainties. 

Figure M shows the excluded regions at the 95% C.L. in the plane (mn±,tan (3) for four values 
of mA, namely 12, 30, 50 and 70 GeV/c 2 , together with the expected exclusion limits. 
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Figure 6: Type I 2HDM: decay branching fractions as functions of the boson masses and tan/3. 
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Figure 7: Type I 2HDM: contours based on the observed p-values CLf, as a function of and tan /3 
or m,A, indicating the statistical significance, N a , of local departures from the background expectation, 
for four values of and two values o/tan/3. The black solid line indicates the change of sign of this 
significance, i.e. where there is a transition from excess to deficit. 
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the charged Higgs boson. The colour-mass correspondence is indicated on the right hand side (units 
are GeY/c 2 ). 
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Figure 9: Type / 2HDM: Excluded regions at 95% C.L. in the (mn±,t&n/3) plane for different values 
of mA-' 12, 30, 50 and 70 GeV/c 2 . The dashed line represents the expected exclusion limit. 
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4 Summary 



The results of the searches carried out by the four LEP experiments for charged Higgs bosons 
predicted by 2HDMs are statistically combined. No significant excess over the SM background 
is observed, and the exclusion limits are extended by several GeV/c 2 with respect to the final 
results of the individual collaborations [2,3,5,6] and the previous combination |T]. In the type 
II 2HDM scenario, assuming that the two decays H + — >cs and H + — >r + u exhaust the H + decay 
width, mass limits are obtained as a function of the branching ratio Br(H + — >t + u). A 95% C.L. 
lower limit on the charged Higgs mass, independent of the boson decay mode, is found to be 
80 GeV/c 2 . A new scenario, for type I 2HDM, thanks to analyses by DELPHI and OPAL in 
the bosonic W*A decay channels, is also performed. In this case, masses of the charged Higgs 
boson below 72.5 GeV/c 2 are excluded at the 95% C.L. for A masses above 12 GeV/c 2 . 
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